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LEAK DETECTION USING 
STRUCTURE-BORNE NOISE 

RELATED APPLICATIONS 

5 

This application claims priority to U.S. Provisional Patent 
Application No. 60/676,428, filed Apr. 29, 2005, hereby 
incorporated by reference in its entirety. 

GRANT REFERENCE io 

The work presented in this application was supported in 
part by a federal grant from the NASA, Grant No. NAG-1- 
02098, the government may have certain rights in this inven- 
tion. 15 

BACKGROUND OF THE INVENTION 

The present invention relates to leak detection and deter- 
mining leak location. To assist in understanding the present 2 o 
invention, problems related to detecting and locating air leaks 
in spacecraft are described. It should be understood, however, 
that the present invention can be used for detection of leaks 
and location of leaks in any number of other types of vessels, 
particularly pressure vessels, and is in no way limited to 2 5 
spacecraft. 

It has been well known for a long time that orbital objects, 
such as meteoroids and space debris, are among the serious, 
but inevitable, threats to spacecraft intended for prolonged 
habitation. Technically, collisions with objects larger than 30 
10 -cm can be avoided using databases compiled by debris- 
tracking systems, such as the Haystack orbital debris radars 
and optical telescopes [1,2]. The strikes of small objects 
between 0.10 and 10.0 mm cannot, however, be avoided using 
the current devices. Furthermore, because the objects strike 35 
the spacecraft at high speeds, up to 10 to 15 km/s [3], even 
small particles in the size range of 1 mm could result in 
penetration of the pressure vessel and subsequent loss of air. 

Referring to several reports [4-7], the major sources of the 
small particles are fragments of payloads, paint flakes, solid 40 
rocket motor slag and discharge, exploded rocket bodies, and 
extraterrestrial micrometeorites. Clearly, there will a higher 
probability of collisions with smaller objects if space explo- 
ration includes trips with longer operational times. 

An example of a current long-endurance spacecraft is the 45 
International Space Station (ISS), now orbiting in low Earth 
orbit (LEO) with an altitude between 3 60 km and 440 km with 
an inclination of 51.6 degrees. It is designed to provide an 
Earth orbiting facility in order to develop advanced technol- 
ogy for human and robotic exploration of space. Ever since 50 
the ISS has launched in 1998, there have been many compo- 
nent replacements, which are mainly caused by the impact of 
high-speed micrometeorites and space resident debris [8]. In 
particular, the impacts of small objects cause air leaks when 
they strike a pressurized module on the ISS. As reported by 55 
NASA in January 2004 [9], parts of the ISS were nearly 
closed off for three weeks due to the air leaks. Although the 
size of the air leak was too small to endanger the crew and 
mission immediately, it was a time-consuming process to find 
and seal the leaks. 60 

Accordingly, immediate location of the source of the leaks 
on the ISS is essential to provide the crew with maximum 
safety, minimize unnecessary mission effort, and maintain 
the operational status. By doing so, the crew can perform 
increasingly ambitious mission goals for longer periods with 65 
less ground supports. However, regardless of size of the air 
leaks, the detection and location of air leaks is a very difficult 
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and time-consuming process because the size of an ISS mod- 
ule and the complexity of its construction. The crew cannot 
inspect the entire series of modules immediately. In addition, 
the ultrasound noise generated by escaping air at supersonic 
speed from small holes into the vacuum of space is not 
audible or detectable inside the ISS because of the nature of 
the leaks. Leaks on spacecraft occur from one atmosphere 
into vacuum and produce almost none of the characteristic 
detectable high-frequency sound, such as that typically 
exploited by industrial leak detection apparatus. 

Conventionally, the crew on the ISS follows the sequential 
module leak isolation process [10]. The process involves very 
inefficient and time-consuming tasks, where the crew moni- 
tors the pressure difference while sequentially closing each 
hatch. Because it is not easy to discern very small pressure 
differences, the isolation process still does not reduce the leak 
detection time and risk of the crew’s safety. Compared to the 
ISS, the MIR, a Russian space station, employs airflow induc- 
tion sensors installed in the hatchways to continually monitor 
airflow and its rate of change [11]. Since the airflow induction 
sensors are designed to detect the very small changes in the 
airflow, they are also sensitive to the circulating air in a 
module. So, the crew must stop all the venting systems in a 
suspicious module to find air-leaks, and this action could 
cause more serious problems. In addition, pinpointing a leak 
location is not possible because the sensors are designed only 
to localize leaking to a specific module losing a measurable 
amount of air, and not to find its exact location within the 
module. 

With the demands of more efficient leak detection systems 
on the ISS, active research has been performed to develop air 
leak detection systems. In 2000, Russian engineers have put 
their best efforts into an acoustic monitor called “OT 2-K” to 
detect air leaks on the MIR [12]. This sensor directs possible 
leak locations by measuring the noise level generated by 
escaping air. In the lab test, the sensor can detect the airborne 
noise in air escaping from a surgical needle at a distance of 
approximately 5 m. However, due to the nature of its working 
principle, it is not an efficient detector or locator of leaks 
inside a module. Instead, the crew must use the sensor outside 
of the module, an extremely time-consuming and often dan- 
gerous activity. Semkin et. al [13] have proposed a multipa- 
rameter transducer, which consists of one ionization sensor, 
four piezoelectric microphones, two thermocouple sensors, 
and a data processing unit in a single package. The dimen- 
sions of the transducer are 80x100x80 mm, and its weight is 
less than 0.35 kg. Analyzing the measurements obtained by 
each sensor, the transducer is capable of detecting gas leaks 
on the MIR from the holes with 0.1-5 mm in diameter at a 
distance of 1 .5 m from the transducer. However, the detect- 
able zone is very small and restricted so that it is not a fully 
functional sensor to locate air leaks on a large-scale space 
station. 

Recently, a company, called CTRL Systems, Inc., has 
developed and tested an ultrasonic leak detector (UL 101®) 
to locate air leaks on the ISS [14]. The detector is a hand-held, 
non-destructive diagnostic tool, which detects ultrasound in a 
narrow frequency band around 40 kHz. The outputs are sup- 
plied to a headset as well as to an analog meter. This method, 
however, is based on the same principle as most of the other 
industrial leak detectors, and is also not efficient for locating 
air leaks on the ISS because it was originally designed to 
detect leaks from pressurized vessels, where air is escaping 
into an atmospheric environment of 1 bar from a pressurized 
vessel at a pressure of at least 1 bar gage, or 2 bar absolute. On 
the contrary, the air leaks on the ISS are leaks into vacuum, 
where the direction of the escaping airflow is into the vacuum 
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of space. Corsaro et al. [15] have proposed a prototype system 
of an acoustic particle impact detector, called “Particle 
Impact Noise-Detection and Ranging On Autonomous Plat- 
form (PINDROP)”. The PINDROP is designed to detect 
hypervelocity impacts by small particles and locate the 
impact sites using the propagation characteristics of the 
acoustic wave generated on a panel. The detector consists of 
a conventional aerogel particle-capture array to collect and 
characterize the small particles in space, acoustic sensor to 
locate their impact sites and autonomous data acquisition unit 
to process the collected data. As a active material in an acous- 
tic sensor, poly-vanylidene fluoride (PVDF) is selected, 
owing to its unique advantages over other piezoelectric mate- 
rials, such as high sensitivity in response to sudden changes in 
strain, very low mass, and flexible installation. To predict 
accurately the signal arrival time, an acoustic propagation 
model is developed and embedded in the data acquisition unit. 
This model is used to locate the impact sites using the relative 
time-of-arrival at three or more sensor locations. Although 
this detector is not designed to find the exact location of air 
leaks on the ISS, it provides the crew with valuable informa- 
tion when they are tracking down suspicious leak locations. 

Kim et. al [11, 16] have developed a leak localization 
software for the ISS. The software uses the measurements 
from pressure gages, spacecraft attitude and rate sensors. 
When air leaks occur, the vent thrust generates a small torque 
on the space station. To preserve the spacecraft’s attitude, a 
reaction torque is needed to stabilize the ISS, depending on 
the size and location of the leak. The software infers from the 
dynamics of the ISS stabilization, considering the geometri- 
cal structure of the ISS, a probable leak location and the 
estimated hole size. According to the developers, their soft- 
ware can locate holes with a diameter as small as 0.4 inches 
(»1 cm). With its outstanding functionalities, NASA is plan- 
ning to employ the software system on ISS in the near future 
[17]. In order to detect and predict the location of air leaks 
accurately, however, exact knowledge of the ISS geometric 
structure is required. If the spacecraft geometry changes, such 
as through the addition of more modules or docking with 
other spacecraft, the software must be entirely reconfigured to 
account for those variations in the spacecraft moment of 
inertia and mass. 

In addition to developing the complex leak detection sys- 
tems, it is desired to utilize the structure-borne noise so that 
the location of a leak is detected by remotely positioned 
sensors. One of the simple and widely used approaches is to 
use cross correlation. The technique assumes in the analysis 
that received signals by sensors at remote locations satisfy the 
conditions of a single mode of propagation with a non-dis- 
persive (or frequency-independent) wave speed. Under the 
assumptions, the received signals are simply non-dispersed 
time-delayed replica of the signal generated at the leak. The 
time delay corresponds to the distance between the sensors 
and the leak source. They are, however, time-shifted with 
different amplitude because of differences in propagation 
paths of the signals. The location of peak on a time scale in the 
cross correlation function corresponds to the time delay, and, 
through the speed of sound in the structure, the location of the 
leak. 

In 1 969, White [18] first used a cross -correlation technique 
to determine a leak location on a thin plate, assuming that the 
structure-borne noise is propagating in all directions with the 
same propagation speed. The temporal shape of the measure- 
ments is also assumed to be unchanging with distance. How- 
ever, the results of the cross correlation do not clearly show a 
large peak. Instead, it shows a gradual rise and fall of sinu- 
soidal function. Ziola et. al. [19] demonstrates a simple leak 
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location method using a very thin aluminum plate, in which 
the method utilizes the classical plate theory and cross-cor- 
relation technique. Given a priori knowledge of the wave 
speed in a certain frequency range, the method can properly 
5 locate the leaks in the thin plate. Their study is limited to very 
thin plates whose thickness is less than 2.5 mm. 

Although a simple cross correlation method is marginally 
capable of locating the leaks, it provides unambiguous indi- 
cations of a leak only if it satisfies the simplified assumption 
to of a single mode of propagation with a non-dispersive (or 
frequency-independent) wave speed. In fact, the structure- 
bone noise in a plate-like structure, such as the outer skin of 
spacecraft, is evidently carried by multiple dispersive modes 
of plate wave propagation. In a dispersive medium, the struc- 
1 5 ture-bome noise propagates with a different velocity at every 
frequency so that the cross correlation technique is not suit- 
able to adequately provide an exact location of the leak 
source. The method provides much poorer results as the 
received signals become more dispersive. Examples of 
20 ambiguous indications in dispersive media are shown in Ref- 
erences 20-23. In addition, the temporal shape of the propa- 
gating structure-borne noise does not remain in constant. 
Furthermore, the effect of the dispersion becomes more seri- 
ous as the distance between the sensors at remote locations 
25 and the leak source increases. As a result, the cross correlation 
technique is not a viable option to locate leaks in dispersive 
media or structures. 

In order to compensate the dispersion-related drawbacks of 
the cross correlation technique in source location, Rewerts et 
30 al. [20-21] have demonstrated a dispersion compensation 
method in a 1 -dimensional structure for water- filled pipe- 
lines. This method isolates particular modes from the dis- 
persed measurements as well as utilizing the advantages of 
the cross correlation technique. Steri et al. [22-23] have 
35 extended the work of Rewert, et al. and showed a robust 
mode-compensating leak location method in a highly disper- 
sive structure, with a leak and sensors arranged in a collinear 
fashion. This method isolates the contributing modes of 
propagation and determines the compensation for the fre- 
40 quency dependence of the contributing mode wave speed. To 
adequately isolate the particular modes, it performs temporal 
and spatial Fourier transforms on the received signals. These 
works also emphasize that the signals need not be collected 
simultaneously for the cross correlations among several pairs 
45 of sensors, owing to the stationary property of the power 
spectrum of the cross -correlated signals. 

Despite these works, significant problems remain in leak 
detection. What is needed is a sensitive and reliable means to 
locate an air leak in pressurized, habitable and long-endur- 
50 ance spacecraft. Such a method should be able to locate air 
leaks accurately with the requisite speed. 

SUMMARY OF THE INVENTION 

55 Therefore, it is a primary object, feature, or advantage of 
the present invention to improve over the state of the art. 

It is a further object, feature, or advantage of the present 
invention to provide for rapid detection of air leaks in a 
spacecraft skin, such as those caused by impact with space 
60 debris or meteorites. 

Another object, feature, or advantage of the present inven- 
tion is the detection of air leaks in a spacecraft skin by using 
ultrasonic waves in the spacecraft skin and not just airborne 
noise. 

65 Yet another object, feature, or advantage of the present 
invention is to provide a continuous system for leak detection 
and the determination of leak location. 
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A further object, feature, or advantage of the present inven- 
tion is to provide for a method and system for detection of air 
leaks that is capable of acconunodating multiple dispersive 
modes of plate wave propagation. 

A further obj ect, feature, or advantage of the present inven- 5 
tion is to provide a framework for leak detection of a pressure 
vessel flexible enough to allow a number of different of sensor 
configurations and corresponding specific methodologies to 
be used. 

A still further object, feature, of the invention is to provide to 
a method for leak detection and location suitable for use with 
any number of types of pressure vessels. 

One or more of these and/or other objects, features, or 
advantages of the present invention will be apparent from the 
specification and claims that follow. 15 

In one embodiment, a coupled phased-array method is 
used. This method provides an all-encompassing data set to 
represent what could be measured if time and complexity 
were irrelevant. Such a method is impractical due the number 
of sensors and complexity of calculations required. However, 20 
use of this method is helpful in fully understanding the prob- 
lems and issued involved. 

In another embodiment, a two-dimensional phased array 
method is provided. This methodology reduces the quantity 
of data required in the coupled array measurement to a pair of 25 
independent, two-dimensional array measurements. For each 
array measurement, one sensor is fixed and the other is 
scanned across an array. This reduces the number of correla- 
tions relative to the coupled phase-array method. 

In another embodiment, a distributed sensor is used. In this 30 
method, leaks are located using only correlations between 
signals from a few isolated point sensors. This method com- 
pares the measured correlations from a few arbitrarily distrib- 
uted sensors with synthetic correlations calculated from all 
possible leak locations and identifies the closest match as the 35 
location of the leak. 

In another embodiment, a dual-sensor method is used. In 
this method, leak detection is performed by measuring leak- 
into- vacuum cross-correlations of noise signals from two 
adj acent transducers as the transducer pair is rotated to deter- 40 
mine the relative phase delay as a function of rotation angle. 
The direction to the leak is found from the variation of phase 
with angle of from synthetic aperture analysis. The leak is 
then located through triangulation from two or more sensor- 
pair locations. 45 

According to one aspect of the invention, a method for 
detection of air leaks in a pressure vessel is provided. The 
method includes sensing structure-borne ultrasound wave- 
forms from a plurality of sensors and determining when cross 
correlation between the waveforms indicates presence of a 50 
leak. The ultrasound waveforms are associated with turbu- 
lence caused by the leak. The method also provides for locat- 
ing the leak using the cross correlation between the wave- 
forms. The ultra sound waveforms may be in a range of about 
300 kHz to 600 kHz. The present invention provides numer- 55 
ous configurations of sensors to be used along with corre- 
sponding methodologies as may be appropriate for different 
environments or applications of the invention. For example, 
the sensors may be configured as a two-dimensional phased 
array, distributed on an inner surface of the pressure vessel, or 60 
may include two rotatable sensors. 

According to another aspect of the invention, a system for 
detection of leaks in a pressure vessel is provided. The system 
includes a plurality of acoustic sensors adapted for sensing 
structure-borne ultrasound waveforms associated with turbu- 65 
lence caused by a leak in the pressure vessel and an intelligent 
control, such as a computer or other computing device, opera- 
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tively connected to the plurality of acoustic sensors. The 
sensors may be of different types, configured in different 
manners and corresponding methodologies are then used for 
leak detection and leak location. The pressure vessel may be 
a spacecraft. 

According to another aspect of the present invention, a 
method for location of an air leak in a pressure vessel includes 
receiving a leak noise signal associated with the air leak in the 
pressure vessel at a plurality of sensor locations to provide a 
plurality of leak noise waveforms, cross correlating the leak 
noise waveforms, and determining location of the leak based 
on results of the step of cross correlating the leak noise wave- 
forms. The results from the cross correlation are a function of 
the leak noise waveform as well as sensor parameters associ- 
ated with the sensors. The sensor parameters can include the 
sensor locations, geometry, etc. A Fourier analysis may be 
applied to the results of the cross correlation to determine 
direction to the leak from the sensors. The location of the air 
leak may then be triangulated based on the direction to the 
leak from each sensor. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram illustrating how a correlation operation 
transforms random noise waveforms into predictable func- 
tion of the leak spectrum, leak location, and detector geom- 
etry. 

FIG. 2 is a schematic of leak and randomly distributed 
sensor arrangement. 

FIG. 3 is a diagram of the mathematical vector notations at 
sensor site A. 

FIG. 4 is a schematic diagram of the measurement configu- 
ration. All dimensions are in millimeters. 

FIG. 5 is a diagram illustrating leak location using mode- 
isolated, 3-dimensional Fourier Transformed cross-corre- 
lated signals at sensor site ‘A’ and ‘B.’ 

FIG. 6 is a schematic diagram of a two-dimensional Piezo 
array sensor system. 

FIG. 7 is a schematic diagram of an 8x8 array with one 
element selected as the reference element. 

FIG. 8 is cross-sectional view of the 64-element piezoelec- 
tric array. 

FIG. 9 illustrates a panel with a leak. 

FIG. 10 is a schematic diagram of a vacuum system. 

FIG. 11 illustrates measured spatial frequency spectra at a 
location A. 

FIG. 12 illustrates measured spatial frequency spectra at a 
location B. 

FIG. 13 illustrates leak location determination by two indi- 
vidual experiments at the test point ‘A’ and ‘B’. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

The present provides for detection of leaks and location of 
the leaks in pressure vessels, such, as but not limited to a 
spacecraft. Although spacecraft are referred to throughout the 
specification, it is to be understood that spacecraft are being 
used merely as one example of a pressure vessel for which 
detection and location determination of leaks is desired. The 
present invention is not to be limited to this specific use. 

The present invention recognizes that, in the case of a space 
craft, listening inside the cabin to locate a leak is ineffective. 
The present invention provides for listening to structure borne 
ultrasonic noise in the skin of the spacecraft. If a leak is 
present, turbulence in the leaking air will generate noise in the 
plate and in the air. Leak noise which couples to the plate-like 
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skin of the spacecraft is carried away from the leak site as 
guided ultrasonic Lamb waves. By monitoring and analyzing 
these waves propagating within the skin of the spacecraft, 
leaks are identified and located. 

The present invention provides for various methods to be 5 
used for leak detection. These methods relay on using mul- 
tiple sensors to sense random structure borne noise wave- 
forms in the skin of a vessel and then applying the mathemati- 
cal operation of cross correlation to compensate for the 
random nature of leak noise. FIG. 1 illustrates a single-fre- to 
quency leak noise signal Ae/ ax propagating in a single plate- 
wave mode to sensors at distance dl and d2 gives a cross- 
correlation that depends primarily on the difference in 
propagation lengths dl-d2. In short, cross correlation trans- 
forms a pair of measured leak noise waveforms into a repro- 15 
ducible function of the leak noise spectrum, geometry, and 
elastic coupling. Even in the more complicated case of mul- 
timode propagation and dispersion (modal velocities that 
vary with frequency), the cross correlation remains a repro- 
ducible function of the leak noise and sensor parameters. 20 

A more detailed mathematical description of various 
embodiments is provided below. These include a coupled 
two-dimensional phased array embodiment, a two-dimen- 
sional phased array embodiment (uncoupled, a distributed 
sensor embodiment, and a two rotatable sensor embodiment. 25 
After disclosure of these particular embodiments, there is a 
separate discussion regarding one of the two-dimensional 
phased array embodiment which is particularly advantageous 
in many circumstances. Then there is a discussion regarding 
a particular sensor system of the present invention. 30 

1 . LEAK DETECTION METHODOLOGIES 

The noise from the leak is inherently random. Moreover, it 
is so faint that it is often buried 20 dB or more under other 35 
noise sources. Therefore, instead of analyzing recorded 
waveforms directly, the statistics of the noise are considered. 

In particular, we want to measure how sound waves detected 
by one sensor relate to sound waves detected by one or more 
additional sensors because this approach yields information 40 
about how the wave propagated to those sensors and hence the 
location of the noise source (the leak itself). One example of 
a preferred statistical tool to do so is the cross -correlation 
function, 

XC ORR (x(t) ,y(t))^r -^(x)y(x-t)dx (1) 45 

The cross correlation at time t of two functions x(t) and y(t) 
is the inner product of x with y delayed by t, as given in Eq. 1 . 
Cross correlation is widely used industrially on leak noise 
generated guided mode waveforms to locate leaks in under- 50 
ground pipelines. Cross correlation transforms a pair of arbi- 
trarily long noise waveforms captured simultaneously at dif- 
ferent sensors into a composite waveform that represents the 
difference in propagation between the leak and each of the 
two sensors. Moreover, the cross-correlation is repeatable 55 
and predictable (ignoring variations from transducer cou- 
pling effects), and recording a full set of all possible cross- 
correlations (including autocorrelations) between sensors 
captures all possible useful information from those sensors 
about the leak. Cross -correlation also compresses extremely 60 
long measured noise waveforms into much shorter wave- 
forms for processing and storage, and transforms very long 
measured noise waveforms with low SNR into shorter high- 
SNR correlations without loss of useful information. 

Consider a single frequency of leak-generated noise in an 65 
infinite plate. Assume the measurable amplitude of the noise 
in guided mode i is I A z -I and its random phase is represented by 
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the complex phase of A z -, which will generally be the same for 
the different modes. The noise in mode i at the location of the 
leak is 


A^xpiflnfl). ( 2 ) 

The waveform that would be measured at a distance d l 
from the leak is 


Aiaitxpijlnft - jk- t di). 


where lq represents the frequency-dependent wavenumber of 
mode i and a z represents combine the distance- and fre- 
quency-dependent attenuative effect of geometric diffraction, 
material absorption, and radiative loss (into the air). The 
cross -correlation between waveforms at distances d a and d 2 is 


^ \M\Mai(di)a t {d 2 )exp( d v - jk t d 2 + jin ft), 
i,i 


or ignoring the i*l cross terms, 


\A\ 2 (*i(di )ai (d 2 )exp [ jk-, (di -d 2 ) + jin ft] 


The correlated noise is no longer random; the cross-corre- 
lation does not depend on the random phase of A z , but only on 
its magnitude. From the point of view of waveform analysis 
the loss factors a z (d x ) and a z (d 2 ) are irrelevant because they 
vary slowly with frequency and hence affect only the broad 
spectrum of the waveform, but not its shape. When we ignore 
the loss terms and cross terms (or assume single mode propa- 
gation) the cross-correlation becomes a spatial function only 
of the difference in path lengths d r d 2 , 


^ \Ai\ 2 expljkiidi -d 2 ) + jlnft\. 


A propagating wave packet that arrives at both sensors simul- 
taneously correlates at t=0, whereas a wave packet that arrives 
at one sensor 1 0 jlls before the other correlates at either t=10jns 
or t=-10 ps depending upon which sensor the wave packet 
passed first. This dependence on the difference in distances 
can be counterintuitive. Because only the difference of dis- 
tances is significant, sound need not come from the same 
location to correlate at the same time. A particular time in the 
cross correlation maps through the modal wavespeed to a 
particular difference of distances. Geometrically, that differ- 
ence of distances occurs along a hyperbole that has its foci at 
the sensor locations, so an observed arrival at that time could 
have come from any location on the hyperbola. 

Successful source location is entirely dependent upon the 
leak noise being localized to a particular region of the plate by 
attenuation, edge absorption, and other losses. In a low-loss 
environment, reflected waves will interfere with direct waves, 
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complicating the process of source location. In the extreme 
case of very low attenuation, resonance or a diffuse sound 
field can occur, and source location will become impossible 
because the cross-correlations will measure nothing but the 
resonance pattern and the Green’s function between the sen- 
sors. Since loss is in general a strong function of frequency, 
only a limited frequency range will contain useful informa- 
tion about the location of the leak. For our tests, we used 
frequencies between 200-600 kHz. Source location in small 
objects requires higher frequency measurement and analysis 
to ensure sound field localization than source location in large 
objects, so lower frequencies may be more useful for actual 
spacecraft than in the smaller test configuration. Of course, 
the present invention is not to be limited to the preferred 
frequency range used for testing. 

Between 200 and 600 kHz and for the 4.76-mm plate 
thickness there are two detectible Lamb modes; the lowest 
order symmetric S 0 compressional mode and the lowest order 
asymmetric A 0 flexural mode. These two modes are both 
dispersive, in that their phase velocity depends on frequency. 
Their dispersion is reflected in the implicit frequency depen- 
dence of the wavenumber k in Eqs. 4-6. The dispersion rela- 
tions k 40 (f) and k^ 0 (f) of the A 0 and S 0 Lamb modes can be 
readily calculated from the known material properties, thick- 
ness and Lamb wave theory. The modes give rise to a total of 
four terms in Eq. 4, two single-mode terms plus two-cross- 
terms. Therefore, the measured correlations will contain dis- 
persed signals from the single mode terms with interference 
from the cross terms. Measured correlations tend to have 
distinct arrivals near t=0 or less distinct dispersed wavetrains 
at t>0 or t<0 depending on the magnitude and sign of the 
difference of distances d 1 -d 2 . 

Cross correlation greatly reduces the amount of informa- 
tion that must be stored. Cross correlations are limited in time 
from the total distances involved and the speed of propaga- 
tion, so they are inherently compact, yet they can represent 
data from arbitrarily long measured waveforms, since the 
integral in Eq. 1 goes out to infinity. Moreover, a full set of all 
possible auto correlations (cross correlation of a waveform 
with itself) and cross correlations between waveforms from a 
set of sensors is complete in the mathematical sense that it 
contains all possible useful information that could be 
extracted from the sensors. Therefore, storing and analyzing 
cross correlations is an effective and efficient alternative to 
storing and analyzing the raw noise waveforms themselves. 

We show by construction that noise waveforms equivalent 
to the originals can be calculated from the full set of auto and 
cross correlations, and therefore that there is no information 
loss from storing and processing only the correlations. Let 
x(t), y(t) be measured waveforms from a pair of sensors, and 
X(f), Y (f) be their Fourier transforms. Since cross correlation 
in the time domain is equivalent to multiplication by the 
complex conjugate in the frequency domain, the cross- and 
auto-correlations are XCORR(x,y)=X(f)Y*(f), XCORR(x, 
x)= X(f)X*(f), and XCORR(y,y)=Y(f)Y*(f). These correla- 
tions can be measured to arbitrarily high signal-to-noise 
ratios by calculating them from sufficiently long recorded 
waveforms. First we reconstruct a frequency domain wave- 
form X r (f) from the known amplitude spectrum VX(f)X*(f) 
and a single random function (|)(f) giving random phase 
The inverse Fourier transform of X r (f) is the reconstructed 
noise waveform x r (t). The reconstructed noise waveform 
Y r (t) is the unique inverse Fourier transform of Y r (f)=[X(f) 
Y*(f)/X r (f)] Since the reconstruction is unique given c|>(f), 
which is equivalent to (but different from) the random phase 
of the leak, the 
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reconstructed waveforms are equivalent to the originals and 
there is no information loss in storing and processing only the 
correlations. 

LI Coupled Two-Dimensional Phased Arrays 

Having now explained the statistical analysis, an embodi- 
ment is now described for coupled two-dimensional phased 
arrays. We have collected an all-encompassing data set to 
represent what could be measured if time and complexity 
were irrelevant. The data set contains all of the possible cross 
correlations between elements of two 2 -dimensional phased 
arrays. Specialization of this data set was used to prototype 
the sparse methods described later herein. 

An automated motion control system was used to position 
15 a pair of 1 .5 -mm diameter piezoelectric transducers used as 
sensors. Computer controlled stepping in 2 mm steps of each 
sensor over its own 32 mmx32 mm grid simulates two phased 
arrays on the 4.76-mm thick aluminum plate. These phased 
arrays have been coupled by measuring correlations between 
20 signals from all possible pairs of locations in the two arrays. 
One-second long waveforms have been recorded at 5 MSPS 
for each pair of sensor locations. Cross-correlations are cal- 
culated from the recorded waveforms and stored. All possible 
cross correlations between the 256 sensor positions in each 
25 array are recorded, leading to a total of 65536 correlations 
arranged as a (1 6x1 6)x(16xl6) phased array calculated from 
655 billion samples. 

The processing of this data to reveal the leak location 
proceeds by straightforward Fourier phased-array analysis, 
30 followed by the application of a-priori knowledge to reduce 
the dimensionality of the data from five dimensions to two. 
The raw correlation data D(t, x 1? y 1? x 2 , y 2 ) is discretely 
indexed along five dimensions: time, x l5 y l5 x 2 , y 2 . A five- 
dimensional discrete Fourier transform converts the data to 
35 frequency -wavenumber space, D(f, k xl , k vl , k x2 , k y2 )- The 
transformed data set is then converted to polar coordinates, 
D(f, IkJ, (|) l5 lk 2 l, c|) 2 ). One dimension can be eliminated from 
the data, along with the cross terms of Eq. 4 by requiring the 
wavenumber magnitudes IkJ and lk 2 l to be equal at the two 
40 sensors, and discarding all data not on the hyperplane 
lk 1 l=lk 2 l. The data can also be reduced in dimension by 
exploiting knowledge of the dispersion relation of the A 0 and 
S 0 modes, k^f) and k so (f), to eliminate data not on the 
hyperplanes Ikl =Uo(f) orlkl=k so (f). The result of this reduc- 
45 tion is a pair of three-dimensional data sets — one for A 0 , one 
for S 0 — each a function of f, 0 l5 0 2 . One more dimension can 
be eliminated by integrating the energy (square of the com- 
plex magnitude) of the data over frequency f, for example 

50 

000(01,92) = f |0(/,fcl,01,fel,02)l l l lH l 2W/(o(/) | 2 ll/ 


55 for the A 0 model This leaves a pair (one for each mode) of two 
dimensional functions of 0 4 and 0 2 that represent the modal 
energy from a common source incident on array 1 at angle 0 : 
and on array 2 at angle 0 2 . 

The A 0 and S 0 mode Eq. 7 energy integrals are then plotted. 
60 The results indicated that the location of the leak is indepen- 
dent of the mode selected for analysis. In an example, the 
estimated leak location is less than 2 mm from the actual 
location. This embodiment, all though exhaustive in scope 
and impractical in terms of time and equipment demonstrates 
65 the capability to estimate the leak location accurately using 
random-signal guided waves. The coupled phased array 
method is particularly insensitive to interference and noise 
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because only a tiny fraction of interference and noise are on 
the extracted hyperplanes; most interference and noise is 
distributed elsewhere in the five-dimensional space and is 
therefore suppressed. 

In summary, the phased array embodiment previously 5 
described uses two or more dense arrays of sensors. Correla- 
tion of one sensor with an entire array is measured. A spatial 
FFT analysis is performed. The method provides a direction 
to the source. The method works in spite of dispersion, multi- 
mode propagation, echoes, etc. However, the amount of data io 
to be collected and analyzed makes this embodiment not as 
practical as other embodiments described later herein. 

1.2 Two-Dimensional Phased Array 

One strategy for reducing the quantity of data required is to 15 
reduce the coupled (x 1? y 1? x 2 , y 2 ) array measurement 
described above to a pair of independent, two-dimensional 
array measurements: (x 1? y x ) and (x 2 , y 2 ). For each array 
measurement, one sensor is fixed and the other is scanned 
across the 16x16 array. The number of correlations required 2 o 
for this measurement is (16xl6)+(16xl6)=512, a factor of 
128 reduction in data from the 65536 correlation full treat- 
ment. The three-dimensional (time, x, y) Fourier transform of 
the data from one sensor array gives a mapping of the direc- 
tions of sound propagation in the vicinity of that array in 25 
terms of frequency and horizontal and vertical wave numbers. 
The three-dimensional transform given measured array 
waveforms D(t,x,y) is: 


D(f, k x , k y ) = J" J* J'j D{ t , x, y)exp(-ik x x - ik y y - i2nft)df dy dx 


To display this as a two dimensional image, we integrate 35 
the magnitude over our selected frequency range 



to obtain the distribution of energy in wavenumber (k x , k^) 
space. 

The bulk of the energy in wavenumber space will be on a 45 
line emanating from the origin. The direction of this line is the 
direction to the leak, and the leak location can be therefore 
found by triangulation from two or more sensor assembly 
locations. The effect of decoupling the two arrays is that 
instead of correlating a specific mode coring one array in one 50 
direction with the same mode crossing the other array in 
another direction, we can see only the angular energy distri- 
bution at each array. No longer is there any guarantee that the 
waves seen at one array match those at the other. In order to 
triangulate we have to assume the same waves are crossing 55 
both arrays, and, for example, in the unlikely event of mul- 
tiple leaks, that assumption might be false. 

Because the required data for the two-dimensional phased 
array measurement is a subset of that required for the coupled 
measurement, instead of performing a separate experiment, 60 
we can reprocess the data from the coupled experiment using 
only the required subset. We arbitrarily select one sensor 
position in sensor array #2 and extract the correlations with all 
possible sensor positions in sensor array #1 to obtain two- 
dimensional phased array data from sensor array #1 D^t, x, 65 
y). Likewise, we select an arbitrary sensor position in array #1 
and extract the correlations with array #2 to obtain phased 
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array data for sensor array #2 D 2 (t,x,y). The data is processed 
using equations 8 and 9. After processing the two dimensional 
phased array data, the rays radiating from the origin indicate 
the measured 0, and 0 2 directions of -8 degrees and 225 
degrees, respectively. Triangulation using these directions 
haves an estimated leak location 2.4 nun from the actual 
position on a 610x610 nun plate. The leak location estimate 
of this method is nearly as accurate as the coupled measure- 
ment. Interference and noise is higher than in the coupled 
case, in part because it has been constructed from far fewer 
data. 

This method has the advantage that it requires far less 
information to find the leak location. It takes less acquisition 
time and less equipment than the full coupled measurement. 
Like the full algorithm, multimode propagation and disper- 
sion will not create interference. Unlike the coupled method, 
this method treats each array independently. It finds the angu- 
lar distribution of wave propagation independently under 
each array and triangulates from the apparently dominant 
direction. This two-dimensional phased array method pro- 
vides a simple, robust method for locating leaks that requires 
far less data to be collected than the coupled phase array. 

1 .3 Distributed Discrete Sensor Method 

While the two-dimensional phased array method requires 
far fewer correlations — and hence far less computation — 
than the four-dimensional method, it still requires two array 
sensors and hundreds of correlations. Using tiny subsets of 
the original coupled data set, we have developed an algorithm 
that locates the leak using only correlations between signals 
from a few isolated point sensors. This algorithm compares 
the measured correlations from a few arbitrarily distributed 
sensors with synthetic correlations calculated from all pos- 
sible leak locations and identifies the closest match as the 
location of the leak. 

The method for locating the leak from the measured cor- 
relations has the following steps: 

1 . Distribute a small number (four in our example) of sensors 
around a plate containing a leak, with known sensor locations 
(Xp yi). ■ • (x 4 ,y 4 ). 

2. Simultaneously record waveforms at all sensors h z (t), 
i=l ... 4. 

3. Measure all (six) possible cross correlations 


Ru(t) = I hi ( r)h t ( t -t) dr 


between noise waveforms from the four sensors, R 12 (t), R i3 
(t) R 14 (t), R 23 (t), R 34 (t). 

4. Calculate dispersion relations k^f) and k so (f) for the two 
fundamental Lamb modes given the known thickness and 
material properties. 

5. Select one mode, A 0 or S 0 . 

6. Select one measured correlation R. 

7. Selection an arbitrary candidate leak location (x c ,y c ). 

8. Calculate a synthetic correlation R s from Eq. 6, ignoring 
the loss factor a, assuming the leak is at the candidate leak 
location, assuming that only the selected mode is present, and 
assuming a flat frequency response over the selected fre- 
quency band. 
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9. Calculate the inner product between the synthetic and 
measured correlations JR(t)R s (t)dt. 

10. Repeat steps 7 through 9 for all candidate leak locations, 
generating a spatial mapping of the magnitude of the inner 5 
product. 

11. Return to step 6 through step 10. multiplying the inner- 
product magnitude for all the measured correlations to create 
a composite mapping of leak intensity for the mode selected 

in step 5 . 10 

12 . Return to step 5 , select another mode and repeat all cal- 

culation through step 1 1 , until all applicable modes have been 
accounted for. Sum the composite mappings to create an 
overall multi-mode map of leak intensity. 15 

The strongest peak in the overall multi -mode intensity map 
provides an estimate of the leak location. 

Because the algorithm processes the different modes inde- 
pendently, it is sensitive to both interference between the 2 q 
modes and interference from cross-terms. For example, a 
correlation signal from an A 0 wave mode but interpreted in the 
S 0 iteration could, when combined with interference in other 
correlations, generate a false peak in the S 0 map. To test the 
robustness of our procedure and to quantify the effect of this 25 
sort of interference, we have run a set of simulations with 
synthetic correlations assuming equal magnitudes of the A 0 
and S 0 modes. These simulations are worst-case in the sense 
that equal A 0 and S 0 magnitudes give the highest possible 
relative amplitude of the cross terms. We performed 32 simu- 30 
lation runs, each with a different, randomly selected, pattern 
of sensors. The average level of the cross -interference peak 
relative to the peak at the leak location is -5 dB. The largest 
interference peak observed is 4 dB stronger than the peak at 
the leak location in the overall map. Under worst-case cir- 35 
cumstances with synthetic data, inter-mode interference can 
generate spurious peaks equal to, or larger than, the peak at 
the leak location. Nevertheless, in each of these cases the leak 
location is also predicted by means of a secondary peak. Since 
the interference is a function of the sensor layout pattern, 40 
careful placement can minimize the interference effects. Add- 
ing sensors beyond the four probes assumed here will also 
tend to reduce interference. Experimental results from several 
sensor configurations with four probes found the actual loca- 
tion of the leak in a 610x610 mm plate with a mean error of 45 
10.4 mm. The arbitrarily distributed sensor method permits 
accurate, rapid leak location with a minimum of pre-posi- 
tioned equipment in an on-orbit spacecraft environment. 

1 .3 Two Rotatable Sensor Method 5Q 

Another sparse-detector method uses two rotatable sen- 
sors. This method is appropriate for cases where a portion of 
the inner surface of the spacecraft pressure vessel is directly 
accessible to astronauts, and is suitable for smaller leaks for 
which the inherent delay involved in a manual measurement 55 
will not be a hazard. By rotating a pair of sensors around a 
central point, correlation data can be collected that can be 
processed by a differential phase analysis or a circular syn- 
thetic aperture analysis to determine the direction of sound 
wave propagation. As in the phased array methods, triangu- 60 
lation from two or more measurement positions will locate 
the leak. 

At each of two or more locations, a pair of rotatable sensors 
is held at a fixed separation and correlations are recorded as 
the sensor assembly is rotated in 15 degree increments. The 65 
direction to the leak is found by examining the phase of the 
Fourier transform of these correlations as a function of angle. 
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Assuming a single mode incident at angle 0 on the sensor pair 
rotated to angle (|), the correlation of the signals from the two 
sensors (separation d) will be 

IAI 2 exp[/'2jt/7 — jkd cos(0 — 0 )]. (10) 

The phase of the Fourier transform of the correlation is 
— j 2 Jtft— jkd cos(0-(|)). The phase varies sinusoidally as the 
sensor assembly is rotated, with the zeroes of the sinusoid 
occurring when the line connecting the sensors is parallel to 
the propagating wavefronts. The extrema of the phase occur 
when that line is perpendicular to the propagating wavefronts. 


P(/) = Yj angle (unwmpped \D{ f, 0))e* 9 ' ‘ ' ’ 

e 

In order for the phase variation to be a meaningful mea- 
sured of leak direction, there must be a single dominant mode 
propagating at each frequency, because multiple modes 
propagating together will cause phase interference. In this 
case, the data indicate that the compressional (S 0 ) mode 
dominates from 450-540 kHz, that the flexural (A 0 ) mode 
dominates from 280-375 kHz, and that other ranges, showing 
large rapid fluctuations, exhibit interference between the 
modes. Once a frequency range has been selected the angle to 
the leak can be estimated by calculating the median phase 
over the selected frequency range of P(f) from Eq. 1 1 , 

medicm [angle{P{f))] (12) 

/ 

Angles calculated with Eq. 1 2 are then triangulated from 
two or more measurement locations to precisely locate the 
leak. 

Another way to analyze these data is as a circular synthetic 
aperture. Unlike the rectangular synthetic aperture (phased 
array) analysis described earlier, the basis functions of the 
circular synthetic aperture are non-orthogonal. Therefore cal- 
culating the angular spectrum from the circular array data is 
an ill-posed inversion problem. To accomplish the inversion 
we use below a modified Lanczos matrix pseudoinverse, 
similar to that discussed by Jackson [14]. 

If a sound wave is incident upon the circular synthetic 
aperture location at angle (|) and the sensor assembly, with 
element separation d, is rotated to direction 0 , then the 
expected correlation (from the simplified correlation of Eq. 6 ) 
would be: 


Y IA | 2 exp[j27r/r - jkidcos(0 - <p )] 


where the sum is over the mode index i. In the synthetic 
aperture analysis we assume incident waves in every possible 
direction and then solve for the amplitudes of those waves 
given the measured correlations. Let A mi be the amplitude and 
phase of the incident wave propagating in the direction in 
mode i. The measured cross correlations with the sensor 
assembly at angle 0 ; would be: 
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XCORRi = \Ami \ 2 exp(j2nt - jkidcosiOt - <p m )) 

m i 


( 14 ) 


5 



df 


( 18 ) 


To solve the synthetic aperture problem, we must solve for 
the angular and modal amplitude spectra IA mi l given the 
measured cross -correlations. Eq. 14 can be represented as a 
matrix E^=[D^ o D, mS J and a vector 10 


2 

1 

2 


15 


such that Eq. 1 4 reduces to the matrix equation 

XCORR^E^C^xpijlnft). ( 15 ) ^ 

Equation 1 5 represents the prediction of correlations from a 
known angular and modal spectrum C q as a matrix multipli- 
cation. Inversion of this equation allows estimation of the 
angular distribution of incident waves from the measured 
data. In a linear synthetic aperture problem, E Iq is a spatial 25 
Fourier transform and is easily inverted. In the current case of 
a circular synthetic aperture, depending on the number of 
chosen values for and 0 l5 this matrix is ill-conditioned and 
may not be square. To construct the pseudoinverse V”. we 
first calculate the singular value decomposition of E, 30 
E=USV r , where U and V are unitary and S is diagonal, with its 
elements S the singular values. Let the largest singular value 
be S mflX . The pseudoinverse is E inv =VS in 'V t , where S*' wv is 
diagonal and constructed from the elements of S, 


5f v 


1 

Sii 


c -°- 2 (^r 


(no sum) 


( 16 ) 


Since E is ill-conditioned, it likely has very small singular 
values. When inverted, these small singular values become 
very large and potentially scale any error or noise in the 
measured correlations. The exponential factor in Eq. 1 6 limits 45 
the noise gain of the inversion process by scaling down the 
inverses of the smallest singular values, following the concept 
of singular value truncation discussed in Jackson [14]. The 
arbitrary factor of 0.2 selects the rate of scaling and corre- 
sponds to the reduction of the inverse of a singular value with 5Q 
half the magnitude of the largest singular value. 

The estimated angular and modal distribution can be cal- 
culated from E lq nv according to 


I 

I AmAQ I 


s C q = ^ EftXCORRte-W 1 

i 


55 

( 17 ) 


Equation 17 gives the synthetic aperture calculation for esti- 
mating the incident angular and modal distribution \A mi \ 2 
from a single frequency component (at frequency f) of a set of 
measured correlations XCORR 1 . By iterating this calculation 
over our frequency range, we obtain the estimated angular 65 
distribution as a function of frequency and mode. Integrating 
this distribution over frequency, 


will yield a peak in the integral in the dominant direction. This 
peak gives an estimate of the direction to the source. 

A two-sensor leak location by tri angulation was performed 
from three measurement positions with both analysis meth- 
ods. The source location error in this case was 2.0 mm for the 
phase comparison method and 2.3 mm for the circular syn- 
thetic aperture method. As before, this source location was 
accomplished in a 610 mm square 4.76 mm thick aluminum 
plate with a 1-mm diameter leak. The two sensor method is a 
viable method for locating leaks. It has the disadvantages of 
requiring manual operation and high sensitivity to cross - 
mode interference, yet it requires a minimum of equipment. 
The equipment need not be manufactured as part of the space- 
craft. 

Each of these methods developed and demonstrated here 
can successfully and repeatably locate air leaks into vacuum 
using structure-borne noise. Of the four methods presented, 
only the first method, the four dimensional phased array is 
impractical because of the extreme amount of data that must 
be collected. Each of the other methods is a viable alternative 
depending on the circumstances involved. A more detailed 
example of the two-dimensional phased array is discussed 
later herein including in conjunction with a miniature array 
transducer with high speed correlator electronics. With a scat- 
tering of two-dimensional arrays on the inside of a spacecraft 
module, this method would instantly identify the location of 
the leak while rejecting possible interference. The distributed 
discrete sensor method would provide the same result, but 
with greater vulnerability to interference and false peaks, yet 
with substantially smaller hardware and computation require- 
ments. The two-sensor method would be appropriate for the 
more spartan spacecraft environment in which the outer pres- 
sure walls are accessible from the inside and for which the 
weight of a permanently embedded monitoring system would 
be prohibitive. 

2. LEAK DETECTION USING 
TWO-DIMENSIONAL PHASED ARRAY 

The concept of two-dimensional phased arrays has already 
been discussed. In this particular embodiment of a two-di- 
mensional phased array, a sophisticated signal -processing 
algorithm is used to determine a leak location. The algorithm 
accounts for the multi-mode dispersive characteristics in a 
plate-like structure, and utilizes structure-borne noise gener- 
ated by turbulence at an air leak. In order to reduce the amount 
of data required and computation efforts, the proposed algo- 
rithm determines only the propagation direction of an air leak, 
rather than time-of-flight information. When an air leak and 
sensors are positioned in a line, the time-of-flight information 
is easily determined by simplifying the computational steps 
in the proposed algorithm. Both cases of determining either 
propagation direction or time-of-flight information are pre- 
sented in this section. We also show the performance of the 
proposed algorithm, using cross -correlated signals on a two- 
dimensional grid at pairs of sensor locations. 

2.1 Mathematical Descriptions of the Signal Processing 
Algorithm 

The proposed signal-processing algorithm is designed to 
isolate particular modes of propagation within the received 
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signals at two sensor sites. To investigate dispersive propaga- 
tion in a plate-like structure, we assume that an air leak source 
and sensors are randomly positioned. For convenience, we 
use the terms “sensor site” and “sensor” interchangeably. For 
example, sensor ‘A’ means that the sensor is located in a 
sensor site ‘A’ . It is assumed that the sensor site is a 2-dimen- 
sional square grid, such as 16x16 or 8x8 array. FIG. 2 illus- 
trates a schematic of a leak and randomly distributed sensor 
arrangement. A leak 50 is shown as well as a first sensor site 
52 and a second sensor site 54. FIG. 3 provides mathematical 
vector notations at the first sensor site. 

First, consider a single frequency w of structure-borne 
noise generated by turbulence at an air leak that has coupled 
into a particular mode of propagation in a plate-like structure. 
Assuming that the sensor is sufficiently far from the leak, the 
geometric dependence of the received voltage signal at the 
sensor A is proportional to r -1/2 . An appropriate mathematical 
expression for the received voltage signal at a sensor ‘A’ for a 
single dispersive mode i is: 

V A (r, (o)=i\T((o) (^exp(/V(P)^) (2.1) 

where N(oo) is the temporal frequency response of the leak, 
k 4 z (oj) is the frequency dependent wave number, and/is the 
distance from the leak to the sensor ‘A’. The term, A z (co), 
represents the amplitude of the coupling of the leak and 
sensor into a particular mode types i (symmetric or anti- 
symmetric), and the sensitivity of the sensor to the displace- 
ments generated by that mode as the structure-borne noise 
passes beneath the sensor. It is assumed that the coupling into 
the mode is uniform with respect to direction. Since the 
structure-borne noise shows multiple dispersive modes in a 
plate-like structure [20-24], the complete mathematical 
model for Equation 2.1 at a single frequency w is the sum of 
n modes of propagation as shown in Equation 2.2 below. 


V\r A , u) = N(oj)(r A y' n Y l 


Similarly, the mathematical expression for the output voltage 
signal at a sensor B is 


V B (r B , (O) = NiuXr 8 ) 112 ^ 5 m (w)exp(j(^(w)r B ) 

m 

5 

where N(co) is the temporal frequency response of the leak, 
kg™(oo) is the frequency dependent wave number, and r 8 is the 
10 distance from the leak to the sensor ‘B\ The term, B m (oo), 
represents the amplitude of the coupling of the leak and 
sensor into a particular mode m, and the sensitivity of the 
sensor to the displacements generated by that mode as the 

structure-borne noise passes beneath the sensor. 

15 

Again, assuming the distance between a leak and sensor is 
sufficiently large, we can consider that a measured wave front 
within a neighborhood of a sensor position at sensor site is 
20 linear. The distance from a leak to sensor ‘A’ can be expressed 
as 

*r A0 + d A -(x A -x A0 ) (2.4) 

25 =r M + d A -s A 

where r 40 is the distance from a reference sensor position x 40 
30 at sensor site ‘A’ to the leak, d 4 is a directional unit vector of 
propagation, x A is the distance from an arbitrary sensor posi- 
tion at the sensor site ‘A’ to the leak, and s A is a vector 
connecting the reference position and the arbitrary sensor 
position at the sensor site ‘A’. Similarly, the distance from a 
35 leak to sensor ‘B’ can be expressed 


+ ■(**-**») (2-5) 

40 

= r 80 + d B -s 8 . 

Substituting Equation 2.4 into Equation 2.2, the output 
voltage signal at sensor ‘A’ is expressed as 


V A (s A , to) = N((o)(r A0 + d A -s A ) 1/2 ^ 4Mexp(j^(w)(/° + d A -s A )) 

n 

= iVM(/°) _1/2 ^ A„Mexp(#"M(r 40 +d A -s A )). 


( 2 . 6 ) 

(2.7) 


Similarly, for sensor ‘B’, 


V B {s B , cl> ) = N{io){r^° + d B ■ 5 B ) 1/2 ^ B m ((L>)Qxp(jkg (c^Xr 80 + d B • s B )) 

n 

= NiuXr 80 )- 112 ^ B m ((x))exp(jk.]j (w)(r B0 + d B -s 8 )) 


(2.8) 

(2.9) 
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In Equation 2.7 and 2.9, it is assumed that reduction of the 
signal amplitude owing to diffraction is negligible. 

Given the expressions for the received signals at sensors ‘A’ 
and 4 B’, we can isolate a particular mode of propagation using 
a spatial window function W(s^) and W(s j5 ). The window 
functions are zero outside the identified neighborhood. For 
example, the simplest window function has a constant non- 
zero value over a square area about the fixed reference point. 
Multiplying the spatial Fourier transformed window function 
by Equation 2.7, 


V B (k B , to) - 


(2.16) 


%)0 B m (u>)W(k B 


V A (k A , to) = J W{s A )V A (s A , w)exp(- jk A s A )ds A 

= J' Wes' 1 ) wM(/°r i/2 y A n (w)<A k 


ds A 


( 2 . 10 ) 


= N(o)(r M ) l,2 y /t„(6j)|J' [W(s A ds A C 4*!S<‘ U| '* 0 ) 


where 25 

J = w(k A -k n A (to)d A ) (111) 

30 

Incorporating Equation 2.1 1 into 2.10, 

35 

V\k\u>)= (2.12) 

AlM(r' ,0 r 1/2 £ A„M W(k A -k’‘ A (to)d A )expU(.k n A Uo)r A0 ). 

n 

40 

When kr 4 =k /, (co) d A , Equation 2.12 can be written as 

45 

V A (k p (to)d A , (o)= (2.13) 

/VM(^°r 1/2 ^ A n {oj)W{k p {oj)d A -k n A (u)d A )expUk n A {u)r A0 ) 

n 

50 

And, if the spatial Fourier transformed window function is 
sufficiently wide, 

^ , (F(co)^,(o)=Ar(co))(^ 0 )- 1 / %(co) expt/FXcoy 40 (2.14) 55 

where 


Assuming k 5 =k g (oo) d B , Equation 2.17 can be written as 

V B {k q {a))d B , to) = (2-17) 

B m (.to)W^(.to)d B -^M^JexpU^Mr 80 ) 


If the Fourier transformed window function is sufficiently 
wide, 

U s (^((o))^- s ,(o)=iV(co(r so )- 1/2 ^ 9 ((o) exp(/^((o)r so ) (2.1 8) 

where 


„ „ (1, m = q 

W((k<’((o)-k$(to))d B )al 

(0 , m*q 


Equation 2.18 represents the mode isolated Fourier trans- 
formed signal measured at sensor ‘B’. 

Since the temporal frequency response of the leak source, 
N(co), is a random signal, we can not directly extract useful 
time-of-flight information by using Equation 2.14 and 2.18 
alone. Instead, assuming mode isolated cross-correlated sig- 
nals are available at sensor ‘A’ and ‘B’ positions, we can take 
advantages of cross-correlation techniques. Given the signals 
at sensor position ‘A’ and ‘B’, the mean mode-isolated, cross- 
correlated signal at sensor A and B positions is expressed as: 

( 2 . 20 ) 

where Orepresents the mean value of the signal in the bracket 
and represents the complex conjugate of a Fourier trans- 
formed signal. Substituting Equation 2. 12 and 2.17 into 2.20, 


W((k p (to)-k n A (to))d A ) 


1, n-p 
0, n ± p 


\C AB (k A ,k B ,a >)) = 


PM . V V A„5^(w)exp(j^(w)r A0 - jk%(to)r B0 ) 

■yr AO r Bi) „ „ 


( 2 . 21 ) 


Equation 2.14 represents the mode isolated Fourier trans- 65 W(k A -k n A (to)d A )W*(k B - k^{to)d B ) 

formed signal measured at sensor ‘A’. Similarly, for a 
received signal at sensor B, 
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where P(co)=<N(<n)N*(co)>, and is the mean power spectra of 
the random, stationary noise, N(co), generated by turbulence 
at an air leak. If k A =kJ 1 (w)d A and k s =k^ m (co)d j5 , Equation 
2.21 can be written as 

( C AB ( k\d A , k^d B , co)) = (2.22) 

P(oj) — = ^ A n (u))B* m (co)exp(jk n A ( (x))r A0 - jk % ( w )r B0 ). 


Equation 2.22 yields a 5 -dimensional dispersion relation as a 
function of k^ w , k Ay , k Bx m , k By m , and co. For example, if the 
modes are considered to be the same in the analysis, the wave 1 5 
numbers are expressed as 

(2-23) 

This leads to 20 


\ C AB (k(co)d A , k(co)d B , co)) = 


(2.24) 




A n (cx))B* m (co)exp(Jk(oj)(r A0 - r B0 )) 


25 


This represents a mode-isolated, temporal and spatial Fourier 
transformed signal at sensor ‘A’ and ‘B\ However, Equation 30 
2.21 is driven under the assumptions that all the data are 
collected simultaneously from the sensor site ‘A’ and ‘B\ In 
practice, this is not feasible using even current state-of-art 
instrumentation. 35 

So, instead of using mathematical analysis shown in Equa- 
tion 2.20-22, we can alternatively perform a 5-dimensional 
Fourier transforms after collecting the mean temporal convo- 
lution between pairs of sensor positions at the sites ‘A’ and 
‘ B ’ . The mean temporal cross-correlation signals between all 40 
sensor positions at sites ‘A’ and ‘B’ is expressed as 

(C AB (s A ,s B , co)) =(V A (s A ,co)V B *(s B ,co)> =P(co) (2.25) 

45 

exp (jk n A (co)(r A0 + d A -s A )~ #£M(r*° + d B ■ s B )) 

50 

Performing a windowed 5 -dimensional Fourier transforma- 
tion, Equation 2.25 is the same as Equation 2.22. When the 
modes are the same (i.e. n=m), 

55 

*/(<o)=V"(<o)=*"( <»). (2-26) 

Then, the leak location is then given by 

60 

G ab (s, p A ,p B ) = J( C AB (k n (co)p A , k n (co)p B )) exp-jb^dco (2-27) 


G ab (s , p A , p B ) = J ( C AB (k n (co)p A , k n (co)p B )) exp -^ (oj)s doj 

= fp(co) -. 1 = Y V A n (u)B* m (co)exp(jk n (co) 

J V r*0r*> ^ ^ 

p A ■ r A0 - jk n (co)p B ■ r B0 )H(k n (co)p A , k n (co)p B , co) 
exp(- jk n (co)s)dco 


and 

^(^(co)^^(co)^,co)=r(^ 4 -A"(co)^ 4 -^) 

W*(k B -k n (a)p B -d B ) (2.29) 

Based on the assumptions described in Equation 2.4, 2.5, 2.15 
and 2.19, Equation 2.28 vanishes if p 4 *d A and p^d^. There- 
fore, the unit direction vectors are 

p A =d A (2.30) 


and 


p*=d*. (2.31) 

Therefore, p A ‘r A0 =R A0 and p B -r B0 =R B0 in Equation 2.28, 
where R A0 and R B0 represent the distance from the leak to 
sensors ‘A’ and ‘B’ respectively. Given these conditions, 
Equation 2.28 can be written as 


G Ab (s , d A ,d B ) = 


(2.32) 


P(co) ' 1 - f A n (co)B*(co)&xp(jk n (co)(R A0 - R m - s))dco 

J r A0 r B0 J 


Equation 2.32 will assume its maximum amplitude when 
s=R A 0 -R B0 . The location of the leak is then determined by 
obtaining the value of s, d A and d B , where G„^(s, d A , d B ) is 
the spectral maximum. Therefore, the leak is located in the 
directions of d A from the sensor site ‘A’ and d B from the sensor 
site ‘B\ 

2.2 Practical Application of the Signal-Processing Algorithm 

In case that a leak and sensor are not positioned in a line, the 
process of locating a leak actually becomes much simpler 
because we seek only the directions of the leak, instead of 
extracting time-of-flight information. Once we find the direc- 
tions, we simply extend the line-of-sight in the directions of 
d A and d B tofind the location of a leak. As a result, the inter- 
section of the lines is the estimated location of the leak. 
Typically, in order to utilize Equation 2.32, we need to per- 
form a five-dimensional Fourier transformation in time and 
space. This process requires high computing power and a long 
computation time. Even though it can be performed with 
exceptional computing equipment, we would have to display 
the results in a five-dimensional domain (k Ax n , k Ay , k Bx m , 
k By m and co) afterwards. Clearly, it is not straightforward to 
locate a leak by examining a five-dimensional data set. 

To overcome such difficulties, we can simplify the process. 
Consider that the mean temporal cross-correlation signals 
between all sensor positions at sites ‘A’ and ‘B’ can be 
expressed as 


where p A and p B represent the unit direction vectors. And sis 65 
the difference in the distance between leak and sensor for sites 
‘A’ and ‘B\ Substituting Equation 2.21 into 2.27, 


(C AB (s A ,s B ,co)> =(V A (s A ,co), V B (s B ,co)) 


(2.33) 
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= P{o ))-= ==X Zj ^«( w )^( w ) ex PWM 

v F F n m 

(r A +d A -s A )-jk%((x))(r B +d B -s B )) 

= PM -^L= C* (/ , w)^ (w)exp(j/c^ M 

V r n 

(r A +d A -s A )) 
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examining the results of Equation 2.35 and 2.36 in the hori- 
zontal (k x ) and vertical (k^) spatial frequency plane. 

2.3 Special Case: Coll inear Configuration 

When a leak hole and sensors happen to occur in a line, we 
cannot directly use the algorithm described in Equation 2.33- 
2.37 to identify the direction of a leak from sensor sites. 
Instead, the time-of-flight information is extracted from the 
received signals. Extracting the time-of-flight information is 
readily performed by letting d A =-d B or d B =-d A in Equation 
2.32 So, 


where 


CV, w) = — !=£ P«i(")exp(-j^M(^° + d B -s B )). 


G ab {s ) -- 


( 2 . 38 ) 


15 


( 2 . 34 ) 


PM , 1 - fA n MP"(w)exp(jk ,2 M(/? 40 - R bo -s))doj 

J r A0 r B0 J 


The C * (s 5 ) represents the received voltage signal for multiple 20 
dispersive modes of propagation at sensor site ‘B\ When a 
three-dimensional Fourier transformation in time and two 
transforms in space (x and y) are performed for only s A , we 
obtain 


( C AB (k A , s 8 , co)) = 


( 2 . 35 ) 


( V A (k A , 0})V B *{k B , u))> 


= P(w)-=C*(/, co) 


2 AnMexp (jk n A ((x))r A0 )W(k A -k n A (aj)d A ) 


30 


As a result, Equation 2.35 shows that the three-dimensional 35 
Fourier transformation isolates the modes of propagation 
only at the sensor site ‘A’ . However, the signal contributions 
from the sensor site ‘B’ still remain as superimposed forms on 
the resulting signal. The maximum spectra magnitude of 
Equation 2.35 in the three-dimensional domain (k^ x ”, k 4y n , 40 
and co) leads to the direction of propagation d A . Similarly, we 
can construct a three-dimensional Fourier transformation in 
time and space for only s B , 


45 


(C AB (s A ,k B , co)> 


( 2 . 36 ) 


( V A (k A , co)V B *{k B , oj)> = P(oj) —== D(s a , co) 

■yfrBO 

^ 5^Mexp(-j/c^ , (w)r B0 )W 3i '(^ -k^(co)d B ) 


50 


where 


55 


D(s a , co) = -Ly A„Mexp(#"M(r 40 +d A -s A )) 

-\ l — * 


( 2 . 37 ) 


60 


The D(s^) represents the received voltage signal for multiple 
dispersive modes of propagation at sensor site ‘A’. The maxi- 
mum magnitude of Equation 2.36 leads to the direction of 
propagation d B . Then, by extending the lines in the directions 65 
of d A and d B , the intersection is found to be the location of the 
leak. In practice, the directions can be readily determined by 


Equation 2.38 yields a peak at s=R^°-R^ 0 with an amplitude 
equal to the integral of Equation 2.38 yields the maximum 
amplitude when s=R" 40 -R £() . The location of the leak is deter- 
mined by obtaining the value of s, where G n AB (s) shows the 
maximum spectra. 

2.4 Simulation 

The proposed algorithm for locating a leak is experimen- 
tally evaluated, using actual data collected from a 584 
mmx609 mmx4.76 mm aluminum plate. In order to simulate 
atmosphere leaking from a habitable space-based environ- 
ment into the vacuum of space, an air leak is introduced by 
drawing air through a 0.71-1.04 nun diameter stepped hole 
using a vacuum pump. As shown in FIG. 4, the measurements 
are collected at a first sensor site 52 (sensor site ‘A’) and a 
second sensor site 54 (sensor site ‘B’). Each sensor site 52, 54 
consists of a 16x16 square grid of measurement points with 2 
nun spacing between the adjacent points, which are carefully 
selected to prevent from mode aliasing [22-23] . At each sen- 
sor site 52, 54, mean cross -correlated signals between a pair 
of transducers, one at the first sensor site 52 and the other at 
the second sensor site 54 are collected on the 16x16 square 
grid. 

Once the signals are acquired, they are band-pass filtered in 
the frequency range of 10-500 kHz. A 1 .3 -nun diameter min- 
iature transducer, called a pinducer, is used to collect the 
signals. Since it is far smaller than any of the sound wave- 
lengths of interest, the pinducer is treated as a point detector. 

In order to utilize the proposed algorithm, we take one 
position at the sensor site ‘A’ (or ‘B’) as a reference point. The 
mean cross correlations with signals at the 16x16 square grid 
points in the sensor site ‘B’ (or ‘A’), against the reference 
point, are Fourier transformed in time and in two spatial 
dimensions (x and y). Then, the energy distribution of the 
measured signals is plotted in the k x -k y plane over the fre- 
quency range of 200-500 kHz. In the k^-k^ plane, the direction 
to the maximum energy distribution is the direction to the 
source leak. Instead of utilizing the energy distribution 
approach, we can alternatively perform a similar data analysis 
using the three-dimensional spectra, one time and two spatial 
dimensions (x and y). In this approach, the resulting three- 
dimensional spectra are plotted in the k^-k^-co plane. Then, a 
plane vertically slicing the k^-k^ plane is rotated with respect 
to the temporal frequency (co) axis until a maximum in the 
spectrum is observed on the slicing plane. The rotation angle 
of the slicing plane showing the spectral maximum is the 
direction to the leak. Although both approaches provide the 
same result, we prefer to use the first approach in the analysis, 
owing to its simplicity. 
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In order to estimate the direction to the leak at the site ‘B’, 
the cross -correlated signals are obtained at all 16x16 points at 
the site ‘A’ while a point at the site ‘B’ is fixed as a reference 
point. Then, its energy distribution is plotted in the in the kjsy 
plane, and an estimated direction to the leak is found to be 
134.96° compared with the actual direction of 1 3 6 . 1 6° . When 
the same data analysis is performed on the cross-correlated 
signals at all 1 6x1 6 points at the site ‘B’ by fixing a reference 
point at the site ‘A’, an estimated direction is found to be 
-5.57°, compared with the actual angle of -6.15°. FIG. 5 
illustrates these plots. 

The variations between the estimations and actual values 
are within the precision of the method and small enough to be 
negligible in practice. Once the estimated directions are 
obtained at sites ‘A’ and ‘B’, we readily find the exact location 
of the leak by drawing lines 58, 60 toward the angles at the 
reference point in the sensor site ‘A’ or ‘B\ Finally, the inter- 
section of the lines at point 62 is the estimated location of the 
leak. The error or distance between the predicted leak loca- 
tion 62 and the actual leak location 64 is found to be 7.8 mm, 
which is acceptable. 

According to the experimental results, it is reasonable to 
assume that the leak location algorithm is effective in locating 
a leak in a plate-like structure, when the leak is from within 
the vessel to a vacuum outside. In addition, the sensitive and 
sophisticated algorithm can pinpoint the location of a leak 
within an acceptable experimental uncertainty. 

Utilizing the 16x16 array measurements, we determine a 
minimal measurement grid size needed to accurately locate a 
leak. The results can be, then, directly implemented in our 
proposed array transducer design. We have chosen the mea- 
surement grid size to be 2x2, 4x4, 8x8, and 16x16 with 2-mm 
point spacing. As the square grid size is increased, the indi- 
cation of a leak source becomes clearer and stronger. For 
example, compared with the results for the 2x2 square grid 
size, ones for the 16x16 grid size show much stronger indi- 
cations. In addition, background noise becomes much smaller 
as the grid size increases from 2x2 to 16x16. 

When a 2x2 square grid size is employed to find a leak 
location, it shows two strong indications of a leak source in 
the k x -k y plane in opposite directions. The directions to the 
strong indications are approximately 9.3° and 1 69°, which are 
assumed to be the estimated direction to the leak source at 
sensor site ‘A’ . We speculate that showing two indications is 
caused by aliasing due to undersampling. Neither of the esti- 
mated directions shows accurate direction to the location of 
the leak, compared with the actual direction to the leak source 
of -6. 15°. Based on the results, using the 2x2 grid size causes 
under-sampling in signal processing so that it provides com- 
pletely wrong indications and higher background noise. The 
high background noise can also increase the error in the 
estimate. Therefore, a square array of 2x2 elements is not 
useful for locating a leak, owing to the seriously reduced 
element count in this case. 

When we employ 4x4 square array elements, the perfor- 
mance improves, and only one strong indication and low 
background noise are observed in the k x -k y plane. The esti- 
mated direction is found to be approximately 0°, compared 
with the actual direction of -6.18°. We still observe a very 
weak indication of the leak source in the opposite direction, 
which is caused by the reflections from the plate edges. How- 
ever, the indication is so weak that it can be negligible in the 
analysis. If we use 8x8 square array elements, the estimated 
direction is found to be approximately -5.38°. This result is 
very similar to the estimated direction of -5.57°, obtained by 
using 16x16 array elements. 
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Based on this analysis, it is reasonable to assume that an 
8x8 square array with 2 nun spacing can still provide an 
acceptable estimate of the location of a leak. In addition, 
computation burden and the amount of data required are 
5 substantially reduced. For our practical device, we therefore 
employ an 8x8 square array of transducer elements in our 
proposed array transducer design. 

3. DESIGN AND FABRICATION OF A 
io PRACTICAL ARRAY SENSOR 

A prototype of a two-dimensional PZT (lead zirconate 
titanate) array sensor has been designed and fabricated to 
detect structure-borne noise generated by downstream turbu- 
15 lence at an air leak at a remote location. The sensor is still 
small enough to be permanently attached to the spacecraft 
and has few enough elements to permit rapid data acquisition. 
After the signals are acquired, they are utilized in our pro- 
posed leak location algorithm to find the direction to the leak 
20 source. For the fabrication of the proposed array sensor, we 
have utilized common PCB (printed circuit board) fabrication 
techniques and standard mirco -electronic components. As a 
result, the overall design can be simply modified, depending 
on the requirements. 

25 The fundamental design philosophy of the array sensor is 
to include all active two-dimensional array elements, 
preamps, multiplexers, integrated electronics and mechanical 
components in a single unit mounted to the circuit board. FIG. 
6 shows a schematic diagram of the PZT array sensor 100. To 
30 achieve satisfactory performance with both a lower compu- 
tational burden and a simplified hardware design, we have 
exploited the experimental results obtained. We choose, 
therefore, an 8x8 square array 102 for our reduction to prac- 
tice of this technique. Such an 8x8 square array 102 with 
35 2-mm spacing is shown in FIG. 7. Each element has a diam- 
eter of 1-mm. We have selected element number 29 as the 
reference element for the correlations. 

The prototype PCB 112 includes all the components, 
except a PZT disk 116 and computer 118. All of the electrical 
40 components are surface mounted type and placed on a 3 M x3 M 
PCB 112, and a PZT disk 116 is chosen as an active compo- 
nent to monitor the structure borne leak noise. The PZT array 
sensor 100 includes a 32 channel mux 104 electrically con- 
nected to a preamp 108 and a second 32 channel mux 106 
45 electrically connected to a second preamp 110. Integrated 
electronics 114 are also provided on the PCB 112. The choice 
of the PZT disk 116 is dictated by the need to have high 
sensitivity and its resonant frequency. We choose a 30 - min in 
diameter, 7 -mm thick PZT disk, which is the thickest PZT 
50 disk with a diameter of 30 mm commercially available on the 
market. To provide electrical contract with the plate during 
the experiments, one side of the PZT disk is metalized with 
Silver. However, there is no metal layer on the other side. 
Instead, the square array pattern shown in FIG. 8 is fabricated 
55 on the backside of the prototype PCB 112, and tightly contact 
with the surface of the PZT disk 16. 

During the data collection, the metalized side of the PZT 
disk 116 is in intimate contact with the plate-like specimen, 
which itself is electrically grounded. When an element is 
60 addressed by a computer 118, the element reports the time- 
varying elastodynamic field in the PZT disk 116 at the par- 
ticular element location. The resulting electrical signal is 
cross -correlated with the signal measured by the reference 
element. FIG. 9 shows a 64-element PZT array sensor proto- 
65 type in cross-sectional view during operation. While data is 
being collected, the PZT disk 116 is held in tight contacted 
with the plate 120, and mineral oil is used as a coupling 
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medium to assure high sound transmission. A leak 122 into a 
vacuum system 124 is shown. After all the data are collected, 
they are processed using the leak detection algorithm. 

A prototype of the array sensor has been tested on a 4.76- 
mm thick, integrally stiffened aluminum panel acquired from 5 
Boeing Aerospace. The panel has four 20-mm high and 4.76- 
mm thick stiffeners across the rear side of the panel. The 
schematic diagram of the panel specimen is shown in FIG. 1 0 . 

To simulate atmosphere leaking from a habitable, space- 
based environment into the vacuum of space, an air leak is to 
introduced by drawing air through a 2-mm diameter hole 150 
using a customized vacuum system 152 . The vacuum system 
152 consists of a vacuum-to-plate adapter 154 with O-ring 
156, 1500-liter buffer volume 158 and a vacuum pump 160 
capable of pumping 6.8 liters per minute of air. FIG. 9 shows 15 
the schematic diagram of the vacuum system 152 . 

In order to efficiently demonstrate its capability of locating 
a leak, we select two test points on the panel as shown in FIG. 

10 . Test point A’ 170 is located at 254 mm straight below the 
leak 150 . But there is a stiffener 174 between the leak 150 and 20 
the test point 170 . Test point ‘B’ 172 is located at 98 mm 
below and 148 mm left of the leak 150 . The coordinates of 
each test point are assumed to be the location of a reference 
element in the prototype, which is near the center of the array. 
During the data acquisition process, the waveform at each 25 
element is averaged 20 times, differentiated, band-pass fil- 
tered, and down-sampled from 5x1 0 6 samples/waveform to 
1 .25x1 0 6 samples/waveform. Then, it is segmented into 1000 
samples before cross-correlation. 

Since all elements are integrated in a single 16-mm square 
unit, including the reference element, we obtain the direction 
to the leak at each test point every experiment. This approach 
reduces the computations needed to estimate the direction of 
the leak source. 

Once the mean cross -correlated signals are acquired at 35 
each element, the signals are Fourier transformed in time and 
in two spatial dimensions (x and y). Then, the energy distri- 
bution of the measured signals over the frequency range of 
200-500 kHz is plotted in the k x -ky plane. In the k^-k^ plane, 
the direction to the maximum energy distribution is the direc- 40 
tion to the source of the noise, presumably the leak. 

FIG. 11 shows that an estimated direction to the leak from 
test point ‘A’ is 89.95°, compared with the actual direction of 
90°. The variations between the estimations and actual values 
are within the precision of the method and small enough to be 45 
negligible in practice. As expected, there is also a weak false 
indication in the opposite direction caused by the reflections 
from the plate edges. Although we have observed weak false 
indications, the prototype sensor and the leak location algo- 
rithm show excellent performance even in the presence of the 50 
stiffener, which would make the operation of a conventional 
single-sensor detector nearly impossible. In our experimental 
configuration, the ratio of the panel’s thickness to the stiffen- 
er’s height is approximately 4.2 (=20-mm/4.7-mm). The 
effect of the ratio on the leak location has been investigated 55 
further. 

Estimated directions to the leak source from test point ‘B’ 
are shown in FIG. 12. Compared with the actual direction of 
56.5°, an estimated direction at test point ‘B’ is found to be 
57.6°. The difference between the estimation and actual 60 
direction to the leak is small enough to be considered negli- 
gible for our purposes. When the prototype sensor is rotated 
clockwise through 90° at test point 4 B’ the direction to the 
leak is found to be 144.9°, compared with the actual angle of 
144.6°. This shows marginally better performance than the 65 
previous experiments in which the prototype is not rotated. 

We speculate that it is caused by the misalignments of the 
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reference point in the prototype with the location of test point 
‘B’. In addition to the misalignment, it also occurs because 
the prototype is not fully rotated 90° in the clockwise direc- 
tion. However, we can minimize such unwanted error by 
taking the measurements at more than two points. 

Once the estimated directions are obtained, we readily find 
the location of the leak by extending lines toward the resulting 
directions at test points ‘A’ and ‘B’. The line is simply drawn 
from the origin to the direction of the maximum in the energy 
distribution in the plot. The intersection of the lines is then the 
estimated location of the leak. FIG. 13 shows the location of 
the leak estimated by extending the lines 170 , 172 in the 
direction of 89.5° and 57.6° at test points ‘A’ and ‘B’ respec- 
tively. The two lines 170 , 172 meet at a point 174 which is 
8.59 mm, directly above the actual leak. 

The described array sensor is effective in locating a leak in 
a plate-like structure with the requisite speed and acceptable 
accuracy, when the leak is from within the vessel to a vacuum 
outside. In addition, the sensitive and sophisticated algorithm 
can pinpoint the location of a leak within an acceptable 
experimental uncertainty. Of course, the present invention 
contemplates numerous variations in the design and construc- 
tion of an array system which can be used in conjunction with 
the methodology of the present invention, including differ- 
ences in array size, frequencies used, and other variations. 

Thus, the present invention provides for leak detection and 
the determination of leak location in a pressure vessel through 
sensing structure-borne ultrasound waveforms associated 
with turbulence caused by the leak. The determination of leak 
30 location uses cross -correlations of the sensed data. In all of 
the leak location methods discussed, cross correlation is 
relied upon to compensate for the random nature of the leak 
noise signal. The cross correlation transforms a pair of mea- 
sured leak noise waveforms into a reproducible function of 
the leak noise spectrum, geometry, and elastic coupling. Even 
in the more complicated case of multimode propagation and 
dispersion (modal velocities that vary with frequency), the 
cross correlation remains a reproducible function of the leak 
noise and the sensor parameters. 

The present invention provides for numerous variations in 
the type of pressure vessel used, number of sensors used, 
placement of sensors, and corresponding methodologies for 
leak detection and location including the manner in which 
cross correlation is used. The present invention is not to be 
limited merely to what is specifically described herein and 
one skilled in the art and having the benefit of this disclosure 
will appreciate the broad spirit and scope of the invention. 

4. REFERENCES 

All references listed herein are hereby incorporated by refer- 
ence in their entireties. 

[1] J. L. Foster, A debris avoidance feasibility study for 
robotic satellites , The Orbital Debris Quarterly News, 8(3), 
2-3, 2004. 

[2] E. Stansbery, J. Foster, and C. Stokely, Haystack orbital 
debris Radar measurements updates , The Orbital Debris 
Quarterly News, 9(1), 3-5, 2005. 

[3] R H. Crisco, Historical small debris collision activities , 
The Orbital Debris Quarterly News, 8(4), 6-7, 2004. 

[4] J. C. Liou, Collision activities in the future orbital debris 
environment , The Orbital Debris Quarterly News, 8(4), 8, 
2004. 

[5] R. P. Bernhard and E. L. Christiansen, Orbital debris as 
detected on exposed spacecraft , The Orbital Debris Quar- 
terly News, 2(4), 3-4, 1997. 



US 7,739,899 B1 


29 

[6] W. P. Schonberg, Protecting spacecraft against meteor- 
oid/orbital debris impact damage: An overview, Space 
Debris, 1, 195-210,2000. 

[7] P. Anz-Meador, ISS Space Shuttles examined for debris 
impacts , The Orbital Debris Quarterly News, 6(2), 3, 2001 . 5 

[8] J. Oberg, Crew finds ‘ Culprit ’ in space station leak, 
www.msnbc.msn.com , Jan. 11, 2003, (Date retrieved: Mar. 
\l ,h , 2005). 

[9] S. Holland, R. Roberts, D. E. Chimenti, and M. Strei, 
Two-sensor ultrasonic spacecraft leak detection using 10 
structure-borne noise , Acoustic Research Letters Online 
(ARLO), 6(2), 63-68, 2005. 

[10] National Aeronautics and Space Administration, Inter- 
national space station complex operations emergency pro- 
cedures: all expedition flights , Mission Operations Direc- 15 
torate Operations Division, Aug. 16, 2000. 

[11] J. W. Kim, J. L. Crassidis, S. R. Vadali, and A. L. Der- 

showitz, ISS leak localization using attitude response , 
AIAA Guidance, Navigation, and Control Conference, 
Montreal, CA, August 2001, AIAA paper #01-4272. 20 

[12] A. Zak, Cosmonauts' flight plan: find the leak on MIR , 
www.space.com, Mar. 31, 2000, (Date retrieved: Mar. 2Y d , 
2005). 

[13] N. D. Semkin, A. N. Zanin, K. E. Voronov, and A. A. 
Kirillov, A device for detecting the location of a gas leak in 25 
a spacecraft module, Instru. Exper. Techn., 46(5), 71 1-716, 

2003. 

[14] G. Studor, Ultrasonic Detectors in Space , www.ctrlsy- 
s.com, May 22, 2002, (Date retrieved: Mar. 2Y d , 2005). 

[15] R. Corsaro, F. Giovane, P. Tsou, J. -C. Liou, D. Buzasi 30 
and B. Gustafson, PINDROP — An acoustic particle impact 
detector , The Orbital Debris Quarterly News, 8(3), 3-4, 

2004. 

[16] J. W. Kim, J. L. Crassidis, S. R. Vadali, and A. L. Der- 35 
showitz, International space station leak localization 
using vent torque estimation, 55th International Astronau- 
tical Conference, Vancouver, BC, IAC-04-A.4.10, October 
2004. 

[17] Anon, “Stopping leaks in space”, Machine Design, 75, 

32, 2003. 

[18] P. H. White, Cross correlation in structural systems: 
dispersion and nondispersion waves, Journal ofthe Acous- 
tical Society of America V. 45 pp. 1118-1128, 1969. 

[19] S. M. Ziola, and M. R. Gorman, Source location in thin 45 
plates using cross -correlation, J. Acou. Soc. Am., 90, 
2551-2556, 1991. 

[20] L. E. Rewerts, R. A. Roberts, and M. A. Clark, Disper- 
sion compensation in acoustic emission pipeline leak loca- 
tion, in Review of Progress in QNDE, 1 6A, 427-434, 1 997. 50 

[21] L. E. Rewerts, R. A. Roberts, and M. A. Clark, The role 
of propagation characteristics in acoustic emission pipe- 
line leak location, in Review of Progress in QNDE, 17A, 
501-508, 1998. 

[22] M. Strei, R. Roberts, S. Holland, and D. E. Chimenti, 55 
Leak location in plates: compensation for multimode dis- 
persion, Review of Progress in QNDE, 23B, 1398-1405, 
2003. 

[23] M. Strei, Acoustic source location in plate-like structure, 

thesis, Iowa State University, 2004. 60 

What is claimed is: 

1. A method for detection of an air leak from a pressure 

vessel, comprising: 


30 

sensing structure-borne ultrasound waveforms from a plu- 
rality of sensors; 

performing a cross correlation of the structure-borne ultra- 
sound waveforms; 

determining from the cross correlation between the struc- 
ture-borne ultrasound waveforms the location of the air 
leak from the pressure vessel. 

2. The method of claim 1 wherein the ultrasound wave- 
forms are associated with turbulence caused by the leak. 

3. The method of claim 1 wherein the ultrasound wave- 
forms are in a range of 300 kHz to 600 kHz. 

4. The method of claim 1 wherein the plurality of sensors 
are configured as a two-dimensional phased array. 

5. The method of claim 1 wherein the plurality of sensors 
are distributed on an inner surface of the pressure vessel. 

6. The method of claim 1 wherein the plurality of sensors 
comprises two rotatable sensors, the method further compris- 
ing rotating the two rotatable sensors relative to the pressure 
vessel through electronic actuation. 

7. The method of claim 1 wherein the pressure vessel is a 
spacecraft and wherein the air leak being into a vacuum. 

8. An apparatus for detection of leaks from a pressure 
vessel, comprising: 

a plurality of acoustic sensor arrays adapted for sensing 
structure-borne ultrasound waveforms associated with 
turbulence caused by a leak in the pressure vessel; 

an intelligent control operatively connected to the plurality 
of acoustic sensor arrays, wherein the intelligent control 
being adapted for receiving a structure-borne leak noise 
waveform from each of the acoustic sensor arrays, cross 
correlating the structure-borne leak noise waveforms 
from the acoustic sensor arrays, and determining loca- 
tion of the air leak in the pressure vessel based on results 
of the step of cross correlating. 

9. The system of claim 8 wherein the plurality of acoustic 
sensors are operatively connected to an inner surface of the 
pressure vessel. 

10. The system of claim 9 wherein the pressure vessel is a 
spacecraft. 

11. A method for location of an air leak from a pressure 
vessel comprising: 

receiving a structure-borne leak noise signal associated 
with the air leak from the pressure vessel at a plurality of 
sensor locations to provide a plurality of structure-borne 
leak noise waveforms; 

cross correlating the structure-borne leak noise wave- 
forms; 

determining location of the air leak in the pressure vessel 
based on results of the step of cross correlating the 
structure-borne leak noise waveforms. 

12. The method of claim 11 further comprising a plurality 
of sensor parameters associated with the plurality of sensor 
locations and wherein the results being a function of the 
structure-borne leak noise waveforms and the sensor param- 
eters. 

13 . The method of claim 11 wherein the sensor parameters 
include the sensor locations. 

14. The method of claim 11 wherein the step of determin- 
ing location of the leak based on the cross correlation com- 
prises applying a Fourier analysis to determine direction to 
the leak from each of the plurality of sensors, and triangulat- 
ing the location of the air leak. 



